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We examine synchronization of identical chaotic systems 
coupled in a drive/response manner. A rigorous criterion is 
presented which, if satisfied, guarantees that synchronization 
to the driving trajectory is linearly stable to perturbations. 
An easy to use approximate criterion for estimating linear sta- 
bility is also presented. One major advantage of these criteria 
is that, for simple systems, many of the calculations needed to 
implement them can be performed analytically. Geometrical 
interpretations of the criterion are discussed, as well as how 
they may be used to investigate synchronization between mu- 
tual coupled systems and the stability of invariant manifolds 
within a dynamical system. Finally, the relationship between 
our criterion and results from control theory are discussed. 
Analytical and numerical results from tests of these criteria 
on four different dynamical systems are presented. 
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If several identical dynamical systems travel the 
same phase space trajectory at the same time 
then they are said to be synchronized. Examples 
of this type of behavior go back at least as far as 
Huygens. Surprisingly, this type of behavior has 
also been observed for chaotic systems. If the 
systems are chaotic then they must be coupled 
together in some fashion. A typical procedure 
for determining a type of coupling that produces 
synchronous motion is to try different couplings 
until one of them works. This paper presents a 
rigorous and an approximate criteria that, if sat- 
isfied, guarantees linearly stable synchronous mo- 
tion for the systems. To use the criterion one only 
needs the equations of motion and the trajectory 
one wants the systems to follow. The paper also 
presents a step by step procedure that will allow 
one to design couplings that satisfy these criteria. 
Since the procedure does not use trial and error, 
and (for simple systems) many of the steps can 
be performed analytically, it is an improvement 
over current practice. 



I. INTRODUCTION 

Since Fujisaka and Yamada's 1983 paper on syn- 
chronized motion in coupled chaotic systems many re- 
searchers have discussed the stability of this type of mo- 
tion Q|. Rather than attempt to reference a complete 



list of papers on this topic we refer the interested reader 
to Refs. §-|. (These papers have extensive bibliogra- 
phies and represent a reasonable introduction to the lit- 
erature.) The discussion in this paper will center around 
the type of synchronization discussed by most of these 
authors. Namely, two or more identical chaotic systems, 
coupled in a drive/response manner, which exhibit mo- 
tion that is chaotic and identical in time. (Although the 
dynamical systems we examine are chaotic our results are 
applicable to nonlinear systems that exhibit other types 
of motion.) 

The central question addressed in this paper is: "Given 
two arbitrary identical dynamical systems how can one 
design a physically available coupling scheme that is 
guaranteed to produce stable synchronized motion." De- 
spite the large amount of effort devoted to this issue there 
are relatively few rigorous results. In most cases rigorous 
results are obtained using Lyapunov functions [^0|-^3[. 
Unfortunately, this method is not regular since, in prac- 
tice, it can only be applied to particular examples. I.e., 
given an arbitrary dynamical system it is not clear how 
one can derive a coupling scheme and/or a Lyapunov 
function which guarantees stable synchronous motion. 

Another rigorous approach is that of Ashwin et. 
al. [|l4j . In a series of papers they presented results that 
are related to the central question. In their work the 
measure (equivalently, the trajectory) used to describe 
the dynamics of the system is of central importance. To 
apply the approach one must show that all normal Lya- 
punov exponents are negative for all measures of the dy- 
namics. Unfortunately, it is well known that Lyapunov 
exponents can be difficult to calculate. Non-rigorous re- 
sults that employ Lyapunov exponents to determine the 
stability of measures to perturbations have also appeared 
in the literature @,|,|l5p0. 

A third rigorous approach by Walker and Mees has re- 
cently appeared ]l8| . They uses the method of Lyapunov 
to develop a sufficient condition for stable synchroniza- 
tion. Their condition relates eigenvalues of the linear 
part of the vector field to the Jacobian of the nonlinear 
part. The analysis in this paper also d differs from that 
in other papers in that they interpret synchronization as 
an observer problem. 

Finally, there are a few special cases where, due to the 
coupling between the nonlinear systems, rigorous analy- 
sis of the stability of synchronization is straightforward. 
One case is when the coupling transforms the driven sys- 
tem into a stable homogeneous linear system with time 
dependent external forcing, (see Ref. |19|). A second is 
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when the coupling between all of the variables is diago- 
nal Jl|,0] . In many practical cases these types of coupling 
can't be achieved. 

The approach advocated in this paper is like that of 
Ashwin et. al. in that it emphasizes the role of indi- 
vidual trajectories within the dynamics of the system. 
However, we do not explicitly use Lyapunov exponents 
or Lyapunov functions to examine the linear stability of 
synchronous motion. It is also like the paper by Walker 
and Mees in that their stability condition and design con- 
cepts are similar the one we develop below. The issue of 
design has also been addressed by Peng et. al. ]l7|] . 

We have two major results. The first is a rigorously 
derived criterion which, if satisfied, guarantees that the 
coupling scheme will yield linearly stable synchronous 
motion on the driving trajectory. The second result is 
a simple, "quick and dirty" criterion that can be used 
to estimate the coupling strength needed for linear sta- 
bility. This second criterion is easy to implement and 
can be used to quickly examine a large range of coupling 
schemes and strengths. For both cases, linear stability 
is with respect to perturbations that are transverse to 
the synchronization manifold. (Of course, because our 
results are obtained from a linearization they can not 
fully address the complications that arise when nonlin- 
ear effects are incorporated. Some of these difficulties are 
discussed in Section M ) 

References g, ||7[|oJ, H and ||, and the results 
presented below, indicate that the stability of the syn- 
chronization manifold depends on the measure used to 
describe the dynamics on the manifold. If the synchro- 
nization manifold contains a set of chaotic trajectories 
then there are an infinite number of possible measures 
for the dynamics on the manifold. Under these circum- 
stances a practical answer to the question of the stabil- 
ity of the synchronization manifold can be quite compli- 
cated. For example, the manifold may be linearly un- 
stable when one uses a measure confined to one periodic 
orbit on the manifold (a Dirac measure full ) and linearly 
stable on a measure confined to a different periodic orbit. 
At the end of the next section we discuss relevant results 
that may allow one to overcome this difficulty. 

We close this section with an outline of the remain- 
der of this paper. In Sections || and III we derive cri- 
teria that indicate when synchronization between drive 
and response systems will occur. The results in Sec- 
tion U are rigorous while those in Section III are ap- 
proximations. Section IV presents results from analytic 
examinations and numerical experiments we have con- 
ducted on the Rossler, Lorenz, and Ott-Sommerer sys- 
tems. In Section ^ we summarize our results and discuss 
how nonlinear effects impact our results. In Appendix ^ 
we discus how our results apply to other systems (in this 
appendix we examine a system of equations that come 
about when considering chaotic masking for private com- 
munications). In Appendix [b] and we present geometric 
interpretations of the manifolds that arise in our discus- 
sion, while in Appendix |c] we discuss a geometrical in- 



terpretation of the theoretical results. In Appendix [l)| we 
discuss the relationship between our results and control 
theory, while Appendix [e| is devoted to an alternative 
analysis of our criterion and an additional examination 
of the Lorenz system. 



II. THEORY PART I: RIGOROUS RESULTS 

This section presents rigorous results regarding linear 
stability of trajectories on the invariant manifold asso- 
ciated with synchronous motion. (Here, and for the re- 
mainder of this paper, when we say synchronization to a 
trajectory is linearly stable we mean it is linearly stable 
with respect to perturbations that are transverse to the 
part of the synchronization manifold defined by the driv- 
ing trajectory.) The major result is a criterion which, if 
satisfied, guarantees linear stability of the driving trajec- 
tory within the synchronization manifold. The criterion 
can be used to design couplings that guarantee linearly 
stable synchronization. 

We will explicitly examine one type of drive/response 
coupling for identical systems that have chaotic uncou- 
pled dynamics p3|-f25t . Assume the dynamics of the driv- 
ing system is given by 



(1) 



where x € M . In the presence of coupling the dynamics 
of the response system becomes 



E(x-y), 



(2) 



where E is a vector function of its argument and repre- 
sents coupling between the systems. We assume E(0) = 
0, hence synchronization occurs on the invariant mani- 
fold given by x = y. Obviously, if the coupling strength 
is below some critical threshold then synchronization will 
not occur. In addition, for some F's and choices of E, 
synchronous motion only occurs within some finite range 
of coupling strengths. For these situations if the coupling 
strength is too small or too large then synchronization 
will not occur. Finally, for some F's and choices of E 
synchronization will never occur. 

Since we are interested in deviations of y from x we 
use Eqs. ([!]) and (Q) to obtain the following linearized 
equation of motion for w = y — x (motion transverse to 
the synchronization manifold) 



dw 
~dt 



[DF(x;t) -DE(0)] w. 



(3) 



In this equation DF(x;t) is the Jacobian of F evaluated 
at x at time t, and DE(0) is the Jacobian of E evaluated 
at w = 0. The synchronization manifold is linearly stable 
if 



lim ||w(t)||= lim ||y(t)-x(i)||=0 

t — >-oo t—*oo 
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for all possible driving trajectories, x(t). (Notice that a 
stability analysis of synchronized chaos needs to consider 
only the invariant trajectories, x(i), that belong to the 
chaotic attractor of the driving system.) We begin de- 
termining the behavior of w(i) in this limit by dividing 
DF(x;t) — DE(O) into a time independent part, A, and 
an explicitly time dependent part, B(x;t), 



DF(x; t) - DE(O) = A + B(x; t). 



(4) 



This decomposition is not unique since we are free to 
add constant terms to A provided we are willing to sub- 
tract the same term from B. This freedom will be resolve 
later when an explicit decomposition is determined. For 
now the reader is asked to accept Eq. (||) as a formal 
decomposition. Denote, and order, the eigenvalues of A 
by 3?[Ai] > 5R[A 2 ] > • • > 5R[A d ], where K[A] is the real 
part of the eigenvalue A. Associated with the eigenval- 
ues are eigenvectors denoted by ei, e 2 , • • • e<j. Next, as- 
sume A can be diagonalized by P = [ei e 2 • • • ej. Thus, 
D = P X AP is a diagonal matrix. Rewriting Eq. M) in 
terms of the coordinate system defined by the eigenvec- 
tors of A yields 



dz 

~dt 



= [D + K(x;t)]z, 



where z = P J w and K = P : BP. Next, define the 
time evolution operator U(t, to) = exp[D(< — to)] and 
use the method of variation of constants to rewrite the 
linearized equation of motion as the following integral 
equation 26 



z(t) =V(t,t )z(t ) + V(t,s)K(s)z(s)ds, (5) 



to 

where K(s) = K[x(s); s] denotes K evaluated along the 
driving trajectory. 

The linear stability of synchronization to the driving 
trajectory, x(t), is determined by the behavior of [|z(t)|| 
in the t — > oo limit. To determine this behavior the 
remainder of the calculations in this section closely follow 
those of Chapter 6 in Ref . |^6| . Begin by using norms to 
obtain the following inequality 1271] 



|z(i)|| < ||U(Mo)||||z(i„) 



||U(M)||||K( s )||[[z( s )||d s . 



It is known that if 3?[Ai] < then ||U(*,t )|| < 
C exp[— fx(t — to)] for appropriate choices of C and fi. In- 
serting this into the inequality and applying Gronwall's 
Theorem pf| yields 



|[z(t)|| <C||z(i )||exp 



(C\\K(s)\\-n)ds 



la 



Hence, a sufficient condition for linearly stable synchro- 
nization is 



fi > lim 



C 



t^oo t — t, 



Jto 



\K(s)\\ds, 



(6) 



which depends on C, fi, ||K||, and the measure of the 
driving signal, x(t). 

The constants C and fj. depend explicitly on the choice 
of matrix norm. A common choice is the Frobenius 
norm lETf 



IMII = 



Tr 



1/2 



a, /3=1 



1/2 



where the f denotes Hcrmitian conjugate and Greek sub- 
scripts denote elements of matrices (M a p is the a f3 ele- 
ment of the matrix M). Applying this norm to U(t, to) 
and using the rank ordering of the A's implies that, in the 
large t limit, ||U(i,t )|| ~ exp[K[Ai](i — to)], an approx- 
imation that becomes exact in the t — > oo limit. There- 
fore, in this limit, C = 1 and /i = — Sft[Ai] . 

Inserting these results into Eq. (^) indicates that the 
sufficient condition for linear stability of synchronization 
along the driving trajectory x(f) is 



U[Ax] > lim 

t— >oo 



1 



t-t 



\K[x( S ); S ]\\ds. 



(7) 



The fact that the right hand side of Eq. (0) is posi- 
tive semi-definite implies that 5R[Ai] must be negative 
for there to be any chance of satisfying the rigorous con- 
dition. 

Because the decomposition given by Eq. ([|) is not 
unique this condition for linear stability of synchroniza- 
tion, Eq. (Q), requires to an examination of K. Notice 

that ||K|| = ||DF - Q|| where DF = P^DFP and 
Q = P _1 [A + DE(0)]P. Thus, the right hand side of 



Eq. (0) can be rewritten as 



I = lim 



t^oo t - t J to 



a 

E (DFaMs);s]-Q af _ 



a, /3=1 



1/2 



ds , 



where the ambiguity of the decomposition is confined to 
the matrix Q. The role of / in Eq. (Q) leads us to con- 
jecture that minimizing / produces the best chance for 
synchronization. We make this conjecture despite the 
fact that 5R[Ai] depends on this decomposition. Hence, 
we define the optimal value of Q = A + DE(O) as the 
one which minimizes /. 

It is straightforward, by setting dl/dQ = 0, to show 
that the optimal value for Q is given by 



Q 



lim 

t— >oo t — to 



DF[x(s);s] ds = (DF) , 



where (•} denotes a time average on the invariant mea- 
sure x(s). This result implies that the optimal decompo- 
sition in Eq. (||) is 
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A= (DF) -DE(O) 
B(x;i) = DF(x;i) - (DF) 



(8) 
(9) 



and the condition for linear stability of the invariant tra- 
jectory in the synchronization manifold is 



5R[A 1 ]>(||p- 1 [B(x;i)]P||) 



(10) 



Equations (|])-([l0]) are the major results of this sec- 
tion. Together they represent definitions and a criterion 
which indicate when synchronous motion along a par- 
ticular driving trajectory is guaranteed to be stable to 
small perturbations in directions transverse to the syn- 
chronization manifold. As indicated by our examples and 
Appendix ^| this criterion can be used to design couplings 
that are guaranteed to result in stable synchronization. 
Unfortunately, the criterion is only sufficient, not neces- 
sary and sufficient. Thus, one can expect, and our nu- 
merical experiments shown, that it is possible for a cou- 
pling scheme to fail this criterion and still produce stable 
synchronization. This is due, in part, to the fact that 
the derivation of Eq. ( |lC)| ) involved inequalities of norms 
which will tend to overestimate the necessary coupling 
strengths. The relationship between these results and 
control theory are discussed below in Appendix [S]. 

As defined above, the decomposition in Eqs. (|J) and (||) 
is optimal in the sense that it minimizes the right hand 
side of Eq. ([h]). We have conjectured that this gives one 
the best chance at satisfying the inequality. To support 
this conjecture we have the following circumstantial ev- 
idence. If we insert Eq. (|^) into a Volterra expansion of 
Eq. (||) then, to second order, the criteria for linear sta- 
bility is 5R[Ai] < (see the next section). For any other 
decomposition this simple approximate stability criteria 
is only correct to first order. If the driving trajectory is 
a fixed point then B(x;i) = and Eq. ( |l0| ) reduces to 
3?[Ai] < 0. This result is what one would expect from 
a simple linear stability. However, it does not occur for 
other decompositions. To address the full optimality is- 
sue one must compare the size of (||P -1 BP||) to the size 
of the first eigenvalue of A. We are unable to make this 
comparison due to the complicated dependence of these 
eigenvalues on the coupling strengths, DE(0). (The idea 
of using a time average to define A is also suggested in 
Ref. @.) 

Because the stability criterion depends explicitly on 
the driving trajectory the measure associated with the 
dynamics is of crucial importance. An early paper by 
Gupte and Amaritkar p2| used unstable periodic orbits 
as driving trajectories for synchronization. They (and 
others) found that, for fixed coupling strength, synchro- 
nization is stable for some driving tra 
stable for other driving trajectories 
determine a coupling strength that results in a stable syn- 
chronization manifold (i.e., synchronous motion is stable 
on all possible driving trajectories) notice that the right 
hand side of Eq. (|l^) is a time average. 

Recently, Hunt and Ott |^8| examined time averages 
of functions for different measures of the dynamics of a 



cctories and un- 
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chaotic system. They found that, for chaotic systems, 
time averages tend to take on their largest values for 
measures confined to the unstable periodic orbits with 
the shortest periods. (This effect was also observed in 
Rcf. H].) They also found that if an unstable peri- 
odic orbit with a somewhat higher period has the largest 
time average then this value is usually only a small in- 
crease compared to orbits with lower periods. Therefore, 
if synchronization is stable to fixed points and periodic 
orbits (with short periods) located in the chaotic attrac- 
tor, one may consider this as an indication that synchro- 
nized chaotic motion is also stable (an assumption also 
proposed in Ref. (8j). This assumption is supported by 
the Hunt and Ott who claim that a nonperiodic orbit 
will yield the maximum time average on a set of measure 
zero in the parameter space associated with the param- 
eters of F. Finally, this paper presents a crude formula 
for determining how high in period one must look in or- 
der to insure that the largest time averages have been 
found. (The issue of obtaining obtain time averages for 
dynamics on chaotic attractors by weighting time aver- 
ages on the unstable periodic orbits has been addressed 
by Cvitanovic p9[ .) 

In general the stability condition given by Eq. ( [To) ) 
requires that one integrate the norm of a matrix along 
the driving trajectory. As such these equations could be 
difficult to use as an explicit test to determine whether 
a given coupling scheme will produce synchronization. 
(Despite this they are still far more efficient than the 
Lyapunov exponent approximations found in most other 
papers.) However, for simple systems, many of the im- 
portant calculations can be performed analytically (see 
below Appendix |E] and Ref. [H for examples). There- 
fore, the criterion of Eqs. (p|)-po|), along with the results 
of Ref. |28| may well allow one to determine coupling 
schemes that insure linear stability of the entire synchro- 
nization manifold. 



III. THEORY PART II: APPROXIMATE 
RESULTS 

Even with the rigorous results of the previous section 
one often desires a "quick and dirty" criterion to approx- 
imate whether a particular coupling scheme and/or cou- 
pling strength will produce linearly stable synchronous 
motion. The major result of this section is just such a 
criterion. 

Begin by noting that an infinite series solution to 
Eq. (||) can be obtained by inserting this equation into 
itself to form a Volterra expansion |30|] 



z(t) 



U(i,io)+^M«(Mo) 



z(i ) 



where 
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M (J) (t, t ) = J dsi J da,- dsj 

x [U(i, Sl )K( Sl )U(s l5 s 2 )K(s 2 ) • • • K( Sj )V( Sj ,to)} . 

Truncating this solution at j = 1 yields the following 
approximate solution to Eq. (jSl) 



(t) ~ [u(Mo)+M«(Mo)l «(to) 



(11) 



(Provided the series is well behaved, it converges, etc., 
one can obtain better solutions than the one discussed 
below by retaining more terms in the series.) To this 
approximation if all elements of U + decay in time 
then limj^o ll z (^)ll = 0, and synchronization is linearly 
stable. 

It is useful to divide the discussion of M^ 1 -* into an 
analysis of diagonal and off diagonal elements. 

1. (Diagonal elements) 
For these elements 



M« =exp[Aa(*-*o)] E P ^ 

t 



la 



DF^(s)ds - (DF^) (t - t ) 



were the optimal decomposition of Eqs. (||) and (||) 
have been used. For this decomposition the term in 
large square brackets vanishes in the large t limit. 
For any other decomposition this term can be ex- 
pected to grow at least linearly with time. Depend- 
ing on the strength of growth this term could offer 
significant competition to the anticipated exponen- 
tial decay from exp [A a (t — to)]. 

2. (Off diagonal elements) 
For these elements 



a/3 



exp (A a * - Apto) 



x / exp [(A/3 - A a )s] K af3 (s) ds. 
Jt 

For the dynamical systems we are examining the 
driving signal, x, is bounded and the Jacobian is 
bounded (typically x is confined to a compact at- 
tractor). Therefore, 



K 
< — 

ap Ap 



(max) 
a/3 



A f 



[exp (Apt) - exp (A a t)] 



This analysis of indicates that if 5R[Ai] < then 
the off diagonal elements of the matrix U + die 
off exponentially for t > to, while the diagonal elements 
die off at a rate controlled by the convergence of the 
time average (DF). Therefore, to second order we have 
the following simple condition for linear stability of a 
trajectory on the synchronization manifold to transverse 
perturbations 



»[AJ < 0. 



(12) 



Equation 
and Eq. 



13) is the major results of this section. It, 
represent a "quick and dirty" criterion to 
determine whether or not synchronization will occur for 
a particular coupling scheme. 

Equation ( |f2| ) is completely consistent with the rigor- 
ous results of the previous section. As expected the ap- 
proximate criterion given by Eqs. (^|) and ( |l2] ) depends 
explicitly on the measure of the driving trajectory. An 
appealing aspect of Eq. ( |f2| ) is that the analytic and com- 
putational burden associated with computing A and it 
eigenvalues is relatively minor. Thus, one can quickly 
search over many types of coupling schemes. 

Our numerical experiments show that Eq. ( |l2| ) is ac- 
tually a good approximation to the minimum coupling 
needed to achieve linearly stable synchronization. One 
possible explanation for the success of this approxima- 
tion is a "folk theorem" from control theory discussed by 
Brogan pl| . "A commonly used stability analysis tech- 
nique is the so-called frozen coefficient method, in which 
all time-varying coefficients are frozen and then the sys- 
tem stability is analyzed as if it were a constant coefficient 
system" M. Bro gan goes on to claim that "when used 
with caution, this approach will usually give the correct 
result" . As a final comment we note that Ref . JTtJ also 
examined A to determine the linear stability of synchro- 
nization. 



IV. NUMERICAL EXPERIMENTS 

This section presents results from numerical experi- 
ments we have performed to test the conditions derived 
in Sections || and [II. The results from the first exam- 
ple are mostly numerical and focus on the approximate 
condition. By the last example almost all of the calcu- 
lations are performed analytically and focus on both the 
rigorous and the approximate condition (see Appendix ^ 
for yet another example). When appropriate we will use 
e( c ) to denote critical coupling strengths where the invari- 
ant trajectory on the synchronization manifold becomes 
linearly stable. 



where K^ ax ^ is the maximum value taken by 
K a p(s) on the driving trajectory and we have as- 
sumed the eigenvalues of A are distinct. 



A. Rossler Example 

The Rossler system is the following set of three coupled 
ODE's 
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FIG. 1. Fixed point, period 1 and period 2 orbits of the 
Rossler systems. 



dx 
lit 

dt 

dz 



x + ay 

b + z(x - c), 



where a = 0.2, b = 0.2 and c = 9. For these param- 
eter settings the dynamics of this system has a chaotic 
attractor, and two unstable fixed points located at x± = 
{-ay±,y±,-y±) where 



y± 



c 

2~a 



1± 1 



Aab 



1/2 



Only the fixed point at x_ is near the attractor. 

To simplify matters we use diagonal coupling, E(x — 
y) = diag(ei, e 2 , e 3 ) • (x — y). For this type of coupling it 
is easy to show that Eq. (||) leads to 



-ex -1 -1 
1 a - e 2 

(z) (x) - c - e 3 



Our numerical experiments examine four different 
measures of the driving dynamics. One is the natural or 
Sinai-Bowen-Ruelle (SBR) measure on the chaotic attrac- 
tor. This arises from following an arbitrary trajectory on 
the chaotic attractor. The others are Dirac measures as- 
sociated with the x_ fixed point, and unstable period 1 
and period 2 orbits. The trajectories associated with the 
fixed point and the periodic orbits are shown in Fig. |l|. 

Because the characteristic equation for A is a cubic 
that can not be easily factored for arbitrary values of the 
e's most of the results were obtained numerically. (We 



FIG. 2. The real parts of the eigenvalues of A as functions 
of e for x-driving. A comes from the Rossler system and the 
driving trajectory is the fixed point X-. The dashed line at 
A = is inserted as a visual aid. 



have calculated these eigenvalues using a symbolic ma- 
nipulator to perform the necessary algebra. The expres- 
sions are not particularly enlightening, and will not be 
presented.) 



1. Fixed point measure 

On this measure B(x;£) = 0, so the rigorous and ap- 
proximated criterion are both 5R[Ai] < (a familiar re- 
sult from stability analysis of fixed points found in text 
books). Figure || shows the real parts of the eigenvalues 
of A as functions of e = e\ when the driving trajectory 
is x_ and we only use the first component of x_ as the 
driving signal (e 2 ~ ej, = 0). This type of driving is typi- 
cally called cc-driving. (A dotted line at A = is inserted 
as a visual aid.) The figure indicates that, initially, Ai 
and A 2 are complex conjugate pairs with positive real 
part. However, as e increases the real parts become neg- 
ative at e(°' = 0.19750.... As e continues to increase 
Ai and A 2 become real at e ~ 1.8, with Ai increasing 
and A 2 decreasing. Ai eventually becomes positive at 
e = 4.9975 . . .. The figure also indicates that A 2 and A3 
merge into complex conjugate pairs but soon split back 
into pure reals at e ~ 9. If the second component of x_ is 
used as the driving signal (y-driving) then a similar fig- 
ure is produced, however Ai remains negative for large e. 

The interesting feature of Fig. || is the prediction of 
synchronous behavior for x-driving when e is in the 
range between about 0.2 and about 5. Although no nu- 
merical experiment can prove the correct values of e", 
our experiments indicates that when x-driving is used 
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a bifurcation from unstable (stable) synchronization to 
a stable (unstable) synchronization appears to occur at 
eW = 0.19750... (4.9975...)- 

The results of all of our numerical experiments on the 
Rossler system appear in Table B. The table indicate 
that for both x and y-driving synchronization to x_ be- 
gins at ~ 0.2. For x-driving synchronization ends 
at e( c ) ~ 4.9, while for y-driving it persists for arbitrar- 
ily large values of e. When z-driving is used Ai and 
A2 are complex conjugate pairs with positive real parts 
throughout the entire range of e. Thus, the rigorous and 
approximate condition predict that synchronization will 
not occur at x_ for this type of driving. The numerical 
experiments appear to verified this result. (We remark 
that the rigorous condition predicts, and we were able to 
verify, a range of values for e which will result in synchro- 
nization about x + for z-driving.) 

2. Periodic orbit measures 

On these measures (as well as the SBR measure) 
B(x;t) 7^ 0, so the rigorous and approximate criterion 
are not the same. We will only consider the approximate 
condition and save a discussion of the rigorous condition 
for later examples. 

The matrix A is a function of time averages over the 
driving trajectory. Table || shows numerically deter- 
mined values for (x), (y), and (z) on each of our mea- 
sures. 

Numerical experiments using the measure associated 
with periodic orbits indicate that the eigenvalues of A 
undergo the same type of splittings and mergings found 
for the fixed point case. More importantly for the pe- 
riod 1 (2) orbit SR[Ai] becomes negative at e ~ 0.04 (0.07) 
for both x and y-driving. In addition, for x-driving and 
the period 1 (2) orbit, SR[Ai] becomes positive at e ~ 4.84 
(4.87). When y-driving is used Ai remains negative for 
large e. The approximate criteria predicts a change in 
stability each time SFt[Ai] changes sign. 

The results of numerical experiments to test these pre- 
dictions are shown in Figs. |^. The figures indicate that 
if a;-driving is used and e = 0.53 then the response sys- 
tem remains a finite distance from the period 1 orbit. 
However, if e = 0.54 then the driven dynamics exponen- 
tially converges onto the period 1 orbit. The figures also 
show that if e = 3.8 then the response system appears to 
diverge while it exponentially converges to the period 1 
orbit for e = 3.7. 

As shown in Table |, numerical tests of these pre- 
dictions indicate that the lower (higher) critical cou- 
pling strength is higher (lower) than that predicted by 
the approximate criterion. Thus, the range of coupling 
strengths for which synchronization occurs is slightly 
smaller than that predicted by the approximate crite- 
rion. This is to be expected since the rigorous criterion 
implies that it is not enough that Sft[Ai] be negative, it 



(a) 




10 -4 I , , , , , , y 1 

200 400 600 800 

Time 

FIG. 3. Results of numerical test for critical coupling 
strengths. The Rossler system synchronized to the period 1 
orbit. The period of the orbit is approximately six units of 
time and x-driving is used, (a) output is every 106/96th of 
the period, (b) output is every 53/96th of the period. 
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must be sufficiently negative to overcome the fluctua- 
tions represented by (||P _1 BP||). However, given the 
"quick and dirty" nature of the approximate criterion we 
believe that it yields reasonable predictions for coupling 
strengths that produce synchronization. 



3. SBR measure 

The results of tests conducted on the SBR measure are 
much more complicated that those obtained on the sim- 
ple measures examined above. Sample plots of |jy — x| 
for various values of e and types of driving are shown in 
Figs. |H These fi gures indicate the bursting phenomenal 
observed by previous researchers when discussing syn- 
chronization. The degree to which the response system 
is "synchronized" to the drive system clearly depends on 
the degree to which one is willing to tolerate bursts. 

For example, when y-driving is used Fig. ^b indicates 
that |jy — x|| ~ 10~ 4 can occur for as long as a thousand 
times around the attractor before climbing to ||y— x|| ~ 1. 
(For e = 0.15 we have observed ||y — x|| < 10~ 4 for as 
many as 4000 times around the attractor.) These exam- 
ples confirm that one must be cautious about discussing 
the condition for linear stability of synchronization when 
the SBR measure is used. The values shown in Ta- 
ble | were selected because ||y — x|| steadily decayed to a 
small value and did not significantly increase for the next 
50,000 time steps. In addition we required that ||y — x| 
exhibit the same behavior for values of e close to and 
on the stable side of the threshold. 

However, because SBR trajectories are dense they will 
eventually come arbitrarily close to any unstable region 
of the attractor. The values of listed in Table Q for 
the SBR measure are outside the range listed for the 
period 1 and period 2 orbits. Therefore, if the SBR orbit 
comes sufficiently close to either of these orbits and/or 
stays close for a sufficiently long time then a burst could 
happen |2(| . The results implied by Table | indicate for 
our numerical experiments the SBR trajectory did not 
come sufficiently close the either the period 1 or period 2 
orbits. 



B. Lorenz Example 

The Lorenz system is the following set of three coupled 
ODE's 



dy 

— = rx — y — xz 
dt y 

dz 

Tt =xy-bz, 

where a = 10, b = 8/3, and r — 60. For these parameter 
values the dynamics of the system has a chaotic attractor 



(a) 




2000 4000 6000 8000 10000 



Time 

FIG. 4. Results of numerical test for critical coupling 
strengths. The Rossler system synchronized to a chaotic tra- 
jectory, (a) x-driving with output every one unit of time, (b) 
y-driving with output every two units of time. 
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and, in practice, Ai could be either A + or A z . 
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FIG. 5. Fixed point, period 1 and period 2 orbits of the 
Lorenz system. 



1. Fixed point measures 

On these measures B(x;t) = 0, and the rigorous and 
approximate criterion are both ?R[Ai] < 0. Because, A z — 
— b — e z < for all positive values of e z we focus our 
attention on A±. Analysis of Eq. (|l3|) indicates that if 
diagonal x-driving is used (e2 = £3 = £4 = 0) then the 
critical coupling strength associated with 3?[Ai] = is 



(ei 



is e 



to 



a(r — 1) = 590. If diagonal y-driving is used 
£3 = £4 = 0) then the critical coupling strength 
= r — 1 = 59. Finally, it is easy to shown that if 



diagonal z-driving is used then £1 — £2 — £3 = £4 = and 
3?[Ai] > 0. Therefore, for z-driving the rigorous criterion 
can not be satisfied. The results of all of our numerical 



experiments on the Lorenz system appear in Table [II 



2. Periodic orbit measures 



and three unstable fixed points, one of which is at the 
origin. 

This example is particularly instructive because if the 
coupling is given by 



DE(0) 



£1 £4 
£3 £2 

£ Z 



then many of the calculations can be performed analyti- 
cally. Furthermore, this type of coupling allows the x or 
y variable to drive both the x and y equations. Ott and 
Ding [^3| have shown that this type of coupling is use- 
ful when not all variables are measurable. Their work, 
and the result presented below, indicate that this type of 
driving may produce or guarantee synchronization when 
purely diagonal coupling does not. 

For this dynamical system it is easy to show that 
Eq. (|) leads to 



-a - £1 
-<*>-<* 

(y) 



a — £4 

-l-e 2 -<*> 
(x) -b - £ 2 



We examine the SBR measure and Dirac measures as- 
sociated with the fixed point at the origin, a period 1, 
and a period 2 orbit. The trajectories associated with 
these measures are shown in Fig. |[ Table || shows nu- 
merically calculated values for (x), (y), and (z), on each 
of these trajectories. The eigenvalues of A are 



A z = -b- e z 

-((7 + 1 



A+ = 



£1 + £2 J 



+ 4-0- 



(1 + ^2)] Z 

£ 4 )(r-(z)-£ 3 )) 1/2 , 



On these measures (as well as the SBR measure) 
B(x;i) 7^ 0. An analysis of the approximate condi- 
tion of Eq. (|l2|), in conjunction with Eq. (|l3|), indicates 
that the approximation to the critical coupling strength 
is £1 = o~[r — (z) — 1] ~ 42 for diagonal x-driving, and 
£2 = r — (z) — 1 ~ 4.2 for diagonal y-driving. Numerical 
experiments on the period 1 orbit are shown in Fig |^. 
They indicate that ~ 17 for x-driving and — 4-2 
for y-driving. Although the error in our approximation 
for x-driving may appear large, comparing it to the value 
eS c ^ = 590 needed for the fixed point at the origin indi- 
cates that it may not be very large. 



(13) 



3. SBR measure 

Results of the tests conducted on the SBR measure 
are much more complicated than those obtained on the 
simple measures examined above. Figures ^ exhibit the 
bursting behavior observed when the trajectory of the re- 
sponse system approaches a region of phase space that is 
unstable to perturbations perpendicular to the synchro- 
nization manifold. Because SBR trajectories are dense 
they will eventually encounter any such region on the at- 
tractor. Thus, although Fig. 0a appears to indicate that 
the system will synchronize when £1 > 18 for x-driving 
we know that the fixed point at the origin is unstable 
for this value of t\. Therefore, if the trajectory comes 
sufficiently close to x = then a burst will occur. 



C. Ott— Sommerer example 

The third dynamical system we examine is the Ott- 
Sommerer [^l| set of four nonautonimous ODE's 
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FIG. 6. Results of numerical test for critical coupling 
strengths. The Lorenz system synchronized to a period 1 
orbit. The period of the orbit is approximately one unit of 
time and output is every 19/50th of the period, (a) a>driving. 
(b) y-driving. 



FIG. 7. Results of numerical test for critical coupling 
strengths. The Lorenz system synchronized to a chaotic tra- 
jectory. The output is every 1 /4 units of time, (a) ^-driving, 
(b) y-driving. 
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dx 

^ = -vv x + Ax (1 - a; 2 ) + y 2 + f sin(^) 

d y 

d,Vy 3 

— = -uv y -2y{x-p)- Aky 



(14) 



where v = 0.05, f = 2.3, w = 3.5, k = 0.0075 and 
p = —1.5. This example is interesting for at least two 
reasons: 

• From the point of view of synchronization, it pos- 
sess a rich structure of invariant manifolds. 

• For the case we will examine, all of the important 
calculations can be performed analytically. 

Originally, Ott and Sommerer examined the two di- 
mensional invariant manifold defined by y = v y — in 
the context of riddled basins. Their results indicate that, 
motion on the manifold is chaotic, however, for our pa- 
rameter values, the manifold is unstable. Thus, typical 
motion for this dynamical systems is in JR 4 , where there 
is one chaotic attractor and no other attracting sets. 

As usual, denote the drive and response systems by, 



and 



E(x-y), 



respectively. 

In principle the driving trajectory, x = [x,v x ,y,v y ] € 
M , however, we will consider driving trajectories re- 
stricted to the invariant manifold of the driving system. 
Under these circumstances, x = [x,v x ,0,0] and DF as- 
sumes a block diagonal form. If we use a block diagonal 
coupling matrix, DE, then the linearized equation of mo- 
tion for w = y x decomposes into motion parallel to, 
and perpendicular to, the invariant manifold of the sys- 
tem, 



til 



DF (i) (x)-DE (1 >(0) 



(15) 



where 



and 



DF (i »(x) = 



DE (x) (0) 





gW(x) 



1±) (J-) 

e l e 4 



An equation similar to Eq. ( pH| ) involving T>¥^\ 
DE"", and g^(x) = 4(1 — 3x 2 ) exists for motion parallel 
to the manifold. This decomposition implies, and our nu- 
merical experiments verify, that coupling strengths exist 
for which the response trajectory collapses onto the in- 
variant manifold of the response system but the response 
system is still not synchronized to the driving system. 
When this occurs, the response system is linearly stable 
to perturbations perpendicular to this manifold but not 
to perturbations parallel to this manifold. 

For the remaining discussion we drop the superscripts 
_L and |j , and leave it to the reader to remember that each 
calculation must be performed in both the perpendicular 
and parallel subspaces. Using Eq. ( |lE| ) it is easy to show 
that 



-ei 1-64 

(g) - £ 3 -(f + £2) 



and 

B(x;t) = 
The eigenvalues of A are 





9(?)-{g) 



A+ = 



e 2 j 



± g [(v + e2-e 1 ) 2 +4(l-e 4 )((g)-e 3 )] 



1/2 



(16) 



and the eigenvectors are N±e± = [(1 — 64), ei + A±], 

where N± are the following normalizations N± — 
1/2 

[(1 - e 4 ) 2 + ki + A±| 2 ] 1 . The eigenvalues of A can be 
real or complex, depending on the values of the e's. Us- 
ing the eigenvectors, e±, and the equation for B, one can 
obtain the following simple expression (in each subspace) 

(IIP-MB^PH) 

= (\g-{g)\) 



' (1 


- £4) 


1/2 


' N 2 + N 2 ' 




-A+| 2 _ 




A+A_ 



Note that this equation for (||P _1 [B(x; t)] P||) di- 
verges at the transition between real and complex A±'s. 
From Eq. ( ^g ) it is easy to see that Sft[Ai] increases as the 
term inside the square root increases from zero. Also, 
if A± are complex then — 5J[Ai] can be made arbitrar- 
ily large by increasing ei and/or £2- These observations 
suggest that the best hope for satisfying the rigorous con- 
dition for linear stability of the synchronization manifold 
lies with choosing e's so that A± are complex with imag- 
inary parts that are not too small. 

It is easy to show that if A± are complex then N + = 
A_ and the rigorous condition for linear stability of syn- 
chronization is 



+ e 1 +e 2 >4(| 5 -( 9 )|)C, 



(17) 



and g^(x) = —2(x —p). 



where 
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c 



-g-^r 

(z/ + e 2 - £l ) 2 +4(l-e 4 )«.g) -e 3 ) 



1/2 



(18) 



(A similar, although more complicated expression for C 
can be obtained when A± are real.) 

The major theoretical results of this example are 
Eqs. (p^)-(]l8|), and the conjecture that the e's should 
be chosen so that A± are complex. Together they rep- 
resent an analytic solution to the rigorous criteria for 
synchronization. 

Equation ( jL7|) is interesting because the right had side 
depends on both the coupling strengths and the driv- 
ing trajectory, while the left hand side only depends 
on the coupling strengths. A useful trick is to set 
C = const. > 0. Under these circumstances the right 
hand side of Eq. ( |l7| ) is constant for a given driving 
trajectory and one can increase the value of ei and e 2 
(maintaining constant C) until the rigorous condition is 
satisfied. (As discussed in Appendix |C] this approach to 
designing couplings that satisfy the rigorous criterion is 
quite general.) 



Because Eqs. (|14|) are nonautonimous they do not have 
fixed points. However, periodic orbits and the SBR mea- 
sure still exist. We examine the SBR and the Dirac mea- 
sures associated with a period 1 and a period 2 orbit 
(see Figs. ||). Numerically calculated values for (g) and 
(\g — (g) |) on each of these measures are shown in Ta- 



ble IV 



We now present examples where the rigorous criterion 
is tested for various types of driving. (Our theoretical 
and numerical results are summarized in Table ^.) Our 
analysis is simple and is designed to shown whether or not 
the rigorous criteria can be satisfied. Thus the numbers 
shown in Table [v| are quite large. As we show in Ap- 
pendix |^, the region of parameter space that results in 
couplings that satisfy the rigorous criteria typically has 
positive volume. Thus, with some additional effort we 
could have found lower values for the coupling strengths 
than those listed in Table 0. 



(a) 
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Period 2 
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(b) 
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1. Diagonal driving 

This type of driving uses all components of x as the 
driving signal, €3 = £4 = 0, and t\ = e 2 = e. For this case 
the usable parameter space is the real line M and rigorous 
results have beenpreviously obtained |l],|7j . For this type 
of driving Eq. (|l^) and Table IV indicate that A± are 
always complex on each of our measures. In addition, 
Eq. (|lj) shows that C is independent of e, and its value 
is fixed by the driving trajectory. Inserting these results 
into Eq. ( |l7| ) leads to the following expression for the 
rigorous condition for synchronization 



e > 



2 (Is -<<?>!) 



-4(5) 



1/2 



(19) 



-2 





X 



FIG. 8. Measures of the Ott-Sommerer model, (a) period 1 
and period 2 orbits, (b) the SBR orbit on the surface of 
section given by cot = mod 2tt. 



12 



Since this condition can always be satisfied linearly sta- 
ble synchronization can always be achieved with diagonal 
driving. 



2. Driving via position 

This type of driving uses only the position variables, x 
and y. The simplest example is when £2 = £3 = £4 = 
and the parameter space is again IR. In order to demon- 
strate that the rigorous criterion can not be satisfied for 
this type of driving we define new parameters u = t\ + v 
and w = 1/C. In terms of the new parameters Eqs. (|ri 
and fllS) are 



uw >4(\g-(g) |) , 
-4(g) = {u-2vf + w 2 . 

On the measures we are considering (g) < so these 
equations define a hyperbola and a circle, respectively. 
The circle constrains the values of the new parameters to 
a one dimensional curve. Thus, the dimension of the us- 
able parameter space is the same for new and old param- 
eters. (C is imaginary for value of t\ s IR corresponding 
to points off the circle.) 

Synchronization to the driving trajectory will be lin- 
early stable if any portion of the circle extends above 
the hyperbola. It is straightforward to show that, on 
the measures we are examining, the hyperbola has w > 
[-4 {g)} 1/2 when u - 2v = [-4 (g)} 1/2 . Thus the hyper- 
bole and the circle never intersect and the rigorous condi- 
tion for synchronization can not be satisfied. (This result 
does not mean that this type of driving will not produce 
synchronization. It only means that our analysis can not 
determine a value of t\ that guarantees synchronization.) 

Another variation of this type of driving uses the po- 
sition variables to drive both the position and the ve- 
locity equation (see Ref |^3| ). For this type of driving 
€2 = £4 = and the parameter space, associated with ei 
and £3, is M 2 . Also, Eq. ( fl6l) and Table IV indicate that 
A± are not complex for all values of £1 and £3. However, 
if Eq. (|lq ) is satisfied then A± are complex. Rewriting 
Eqs. (|17|) and Eq. (18) leads to the following expression 
for the rigorous condition for linearly stable synchroniza- 
tion 



ei >-v + 4C{\g-(g)\), 



£3 



1 < \ 2 



(9) + ^ 



(20) 
(21) 



where C > is arbitrary. 

These equations define a line and a parabola, respec- 
tively. The parabola exist for all possible values of e l5 
and all driving trajectories, and its minimum value is 
determined by the arbitrary constant C. The line also 
exists for all driving trajectories. (Of course one should 
choose C so that this line is in the positive t\ half 
plane.) These curves are guaranteed to intersect, and 



synchronization is guaranteed to be linearly stable for 
any point on the parabola whose t\ value is larger than 
the one associated with this intersection. We remark 
that C — 1/4 (\g — (g) |) is a choice which greatly sim- 
plifies the guaranteed synchronization condition (this is 
the form used to generate the numerical values shown in 
Table |v|) 



ei > 1 



£3 



{v-e 1 f+ ( g )+A(\g-(g)\y 



3. Driving via velocity 

This type of driving uses only the velocity variables, 
v x and v y . As an example, let the velocity variables 
drive both the position and the velocity equations. Thus, 
€\ =£3 = and the parameter space is M 2 . Again we 
notice that if Eq. (|l^) is satisfied then A± are complex. 
However, even if Eq. ([l8]) is satisfied, values of £2 and £4 
which satisfy the rigorous condition for synchronization 
do not exit. To see this define new parameters u = ^ + £2 
and w = (£4 — 1)/C In terms of the new parameters, 
Eqs. @ and @ are 

u>4C(| fl -<fl>|) 
(2C(g)) 2 = u 2 + (2C (g) ~ w f , 

where C is arbitrary. These equations define a line and a 
circle, respectively. It is straightforward to show that the 
circle does not intersect the line on the measures we have 
examined. Therefore, the rigorous condition for synchro- 
nization can not be satisfied by this type of driving. 



V. SUMMARY 

In this paper we investigated the linear stability of the 
invariant manifold associated with synchronous behavior 
between coupled dynamical systems. (See Appendix [b| 
for a discussion of these manifolds.) Although our for- 
malism examined a particular type of coupling our results 
are valid for more general types of coupling, and they 
can be used to determine the linear stability of invariant 
manifolds within a dynamical system (see Appendix ^). 

We have two major results. The first is a rigorous 
criterion, given by Eqs. (|^)-(|To|). If it is satisfied then 
linear stability of synchronization to the driving trajec- 
tory is guaranteed. The condition is based on norms of 
deviations from synchronous behavior. As such it tends 
to over estimate the coupling strength needed to achieve 
synchronization. (An alternative formulation of the rig- 
orous criteria can be found in Appendix |e|) . The second 
result, given by Eq. (|l2|), is a "quick and dirty" crite- 
rion for estimating the coupling strengths need to pro- 
duce synchronous behavior. The approximation is easy 
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to calculate and can be used to quickly investigate many 
coupling schemes and strengths. A major advantage to 
our approach is that both the rigorous and approximate 
criterion have a geometric interpretation that can be used 
to design couplings schemes and strengths that will result 
in synchronization (see Appendix Q). Both criterion can 
also be related to fundamental issues of control theory 
(see Appendix . 

Both the rigorous and approximate condition are de- 
pendent on the measure used for the driving dynamics 
and can yield different stability results for different driv- 
ing trajectories. Nonetheless, based on our numerical 
results, and previous work by others, it may be possible 
to use our criterion to determine couplings that result in 
linearly stable synchronization for arbitrary driving tra- 
jectories p8[ . 

To test these criteria we have performed numerical 
experiments on four different dynamical systems, the 
Rossler, Lorenz (see also Appendix |e|), Ott-Sommerer 
system and a system used to model chaotic masking in 
communication (see Appendix |X|). When taken together, 
these examples provide a thorough examination of the 
stability criterion we propose. 

We close with a discussion of how noise and nonlin- 
ear effects influence the conclusions one can draw from 
the linear stability analysis. In this context consider a 
chaotic attractor that contains an unstable fixed point at 
x = 0. Furthermore, assume that x = is the driving 
trajectory, and the response system synchronizes onto 
this fixed point as the coupling strength increases. The 
change in stability as e increases is a bifurcation which 
we will model by a pitchfork, (see Fig [|). A linear sta- 
bility analysis does not take into account the existence 
of the unstable trajectories for e > e' c '. For arbitrar- 
ily small noise amplitude there exists a range of e values 
near where noise will eventually push the response 
system above or below the dashed lines. When this oc- 
curs the response system can not collapse back to the 
synchronized state and is force to seek out an attracting 
state away from the synchronization manifold. Also, a 
linear stability analysis does not take into account the 
existence of the unstable trajectories near x = when 
e ~ e*. If e ~ e* then noise may cause a loss of synchro- 
nization. Therefore, although a linear stability analysis 
may seem to guarantee stable synchronous motion, noise 
and nonlinear effects may prevent long term synchronous 
behavior. 
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APPENDIX A: RELATED PROBLEMS 

With minor modifications the techniques and results 
discussed above are valid for more general types of cou- 
pling. They can also be used to study the stability of 
invariant manifolds of a dynamical system. In part, this 
generality exits because the linearized stability equations 
we examine arise in a variety of other problems, many of 
which have recently appeared in the literature. The key 
point to realize is that a linearized equation, similar to 
Eq. ([|) , arises whenever a dynamical system has dynam- 
ics on a smooth invariant manifold of its full phase space, 
and one is considering the linear stability of this manifold 
to perturbations which are transverse to the manifold. 

For example, the equations examined by Ott- 
Sommerer are of the form 

!=M( X) y;*)y, 

where M is an n x n matrix whose elements are functions 
of x G M m and y £ M n and time, t. For this system 
y = is an invariant manifold. The linearized equations 
of motion for w = y — are the same as Eq. ( p5| ) with 
DE (±) = (see Eqs. (7) and (8) in Ref @). If this 
manifold is unstable then one can use our criterion to 
determine a DE^ that stabilizes this manifold. By ab- 
sorbing DE' 1 ' into the matrix M(x, y; t) one can deter- 
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mine values for the parameters in M that yield a linearly 
stable invariant manifold. 

Another type of problem that is covered by our formal- 
ism is synchronization between mutually coupled identi- 
cal systems. The equations of motion for this problem 
are often of the form 



rfx 
~dt 
dy 
dt 



F(x)+Ei(y-x) 
F(y)+E 2 (x-y), 



where Ei and E2 are functions of their argument and 
Ei(0) = E2(0) = 0. Because of the lack of constant 
terms in Ei and E 2 the synchronization manifold x = y 
is invariant. The linear stability of this manifold is deter- 
mined by Eq. |j) where DE(0) = DEi (0) + DE 2 (0). Us- 
ing our criterion to determine an E that results in linearly 
stable synchronization is equivalent to determining forms 
for Ei and E 2 . (Here, the nonuniqness of the decomposi- 
tion DE = DEi +DE 2 can probably be reconciled using 
details of the specific problem under investigation.) 

Finally, synchronization in drive response systems of 
the form 



dx. 
~dt 

dt 



F(x) 

G(y,x;e) 



where F = G when x = y and/or e = 1 has been pre- 
viously examined p^ , p4p3] . For this system x = y is 
an invariant manifold for all values of e. (We present an 
example of this type of system below.) 

For each of these problems the important research is- 
sues centers around the stability of an invariant manifold 
of the dynamics. In particular, synchronization between 
N coupled dynamical systems can often be thought of 
as motion on a smooth invariant manifold which lives in 
the full phase space of a single large dynamical system 
consisting of the N smaller systems. Given this interpre- 
tation, our results indicate that the stability of synchro- 
nization is a special case of the stability of invariant man- 
ifolds of a dynamical system. (A similar observation can 
be found in Ref. Therefore, although our results are 

presented in the context of synchronization between iden- 
tical chaotic systems coupled in a drive/response manner 
they are applicable to areas of research where the stabil- 
ity of invariant manifolds is the central issue. Examples 
of this are riddled basins and on-off intermittency, syn- 
chronization of identical systems with mutual coupling, 
and generalized synchronization ]14|,pl|,p2|-p6| . 



1. Communication example 

This example is a modification of a synchronization 
method which has previously been used to experimen- 
tally demonstrate communication via modulated chaotic 




FIG. 10. A bl ock diagra m o f electronic circuits that are 
modeled by Eqs. (Al) and (A4). 



signals |32|| . The modification involves introducing a pa- 
rameter, e, whose variation changes the stability prop- 
erties of the synchronization manifold. Introducing this 
additional parameter, and carefully selecting its value, 
can enhance the stability of the synchronization mani- 



fold [ 19| , p3| . Therefore, adding this parameter is very 
useful for applications. 

A block diagram of the modified drive and response 
circuits is shown in Fig. [l^. Each circuit consist of a 
nonlinear converter, N, which transforms input voltage, 
u, into output, af(u) (see Ref. |l^] for details). The pa- 
rameter a characterizes the gain of N around x = 0. The 
nonlinear amplifier has linear feedback which contains a 
series connection to a low-pass filter, RC , and a resonant 
circuit LC. 

Our analysis requires knowledge of the equations of 
motion for the systems. The dynamics of the chaotic 
driving circuit can be described by the following system 
of three ODE's 



dxi 
~~dt 
dyi 
dt 
dz\ 

~dt 



Hi 



-Xi 



5yi 
7[a/(ati + 



-'1 



(Al) 



zi] - cry 1, 



where the parameters a, 7, 5, e and a are all positive. 
(The relations between the parameters of the model and 
the parameters of the circuits can be found in Ref. [Q ) . 
In these equations m denotes the time dependent mes- 
sage one wants to send, although X\ + m is the actual 
transmitted signal. 

In numerical simulations we use 



f(x) = -sign(x) 




-.1/2N 



(ji(x)-a) 2 + C 



(A2) 



a)]/{b- 



if |x| < a 

if a < \x\ < b 

if \x\>b. 



(A3) 
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If we set a = 0.5, b = 1 .8 and c = 0.03 then f(x), as given 



by Eqs. ( A2) and (A3), fits the nonlinearity of the con- 
verter, N, to within w 2% accuracy. The validity of the 



model (Al),(A2) has been confirmed via synchronization 
between the real circuit and the model S. 

The response system is driven by the voltage x\ from 
the drive system. Experimentally, it is fed into another 
circuit which couples it to the transmitted signal, x\ +m 
(see Fig. [!(]). The dynamics of the response circuit is 



dx 2 
~~dt 

dm 

dt 
dz 2 
~dt 



-x 2 



Z-2 



5y2 

= j[af[x 2 + e(xi + m 



(A4) 



*2\ 



The parameter, e, indicates the strength of the cou- 
pling between the drive and response circuit. Note that if 
e = then the drive and response system are uncoupled. 
If e = 1 then a driving signal, x\ + m, is the argument of 
the nonlinearity of the response circuit. If the drive and 
response circuit synchronize then the response system op- 
erates as a "chaos filter" which can be used to extract the 
message from the transmitted signal. (To see this, notice 
that if e = 1, and the response system synchronizes to 
the drive system, then x 2 — x\ and the message, m, can 
be recovered by subtracting x 2 from transmitted signal.) 
More importantly for applications, if the the two systems 
synchronize for e ^ 1 then this same procedure can still 
be used to transmit messages. 



An examination of Eqs. (A.1) and (A4) results in the 
following equations for A and B 



0-10 
-1 -6 1 
—a 7 



and 



B(xi 





ff(xi) 



where <?(xi) = cry(l — e)f'(xi+m) and f'(x) 
The characteristic equation for A is 



df(x)/dx. 



A 3 + (S + 7 )A 2 + (1 + a + 7 S)A + 7 = 0. 

This equation is not easily solved for the eigenval- 
ues. However, the Routh Hurwitz criteria indicates that 
5ft [A] < for all eigenvalues [0. 

For this example, A is independent of e and 
(HP -1 [B(x)] P||) = when e = 1. Therefore, if e = 1 
then the rigorous condition for linear stability of the syn- 
chronization manifold (3?[Ai] < 0) is satisfied. If e ^ 1 
then (||P _1 [B(x)] P||) ^ and this simple analysis fails. 

To understand the importance of e ^ 1 for ap plic ations 
notice that if f(x) is given by Eqs. rtA2) and (A3) then 




FIG. 11. The transverse Lyapunov exponents as a function 
of e. 



f'(x) is not continuous. However, this function is only 
an approximation to the real f{x). If the real f(x) has a 
continuous derivative (in fact we believe f{x) is smooth) 
then there exists an open neighborhood of e = 1 where 
the rigorous condition for synchronization (Eq. ([l0])) is 
guaranteed to be satisfied. 

The existence of this neighborhood is important be- 
cause, although the synchronization manifold is linearly 
stable for e = 1 this does not guarantee that e = 1 re- 
sults in the most stable manifold. Because the loss of 
synchronization is unfavorable for applications involving 
communications one wants a synchronization manifold 
that is as stable as possible, and which recaptures the 
response trajectory as soon as possible in the event that 
synchronization is lost. Therefore, one should to use an 
e value that yields the most stable synchronization man- 
ifold. 

In Fig. [H] we show the transverse Lyapunov exponents 
as a function of e for our model. The figure indicates 
that, at least as far as transverse Lyapunov exponents 
are concerned, the most stable synchronization manifold 
occurs when e ^ 1. The fact that the rigorous condition 
is guaranteed to be satisfied in some neighborhood of 
e = 1 means that we are free to search for values of e 
that result in the most stable synchronization manifold. 



APPENDIX B: GEOMETRICAL DISCUSSION OF 
INVARIANT MANIFOLDS 

In this appendix we discusses the various invariant 
manifolds that arise when discussing synchronization. 
The idea of an invariant manifold is fundamental to our 
approach. Because we have discussed many different in- 
variant manifolds it is useful to carefully construct a ge- 



16 



ometric picture of these manifolds. This construction 
provides a unified interpretation of the issues discussed 
above. Let x 6 M d denote the trajectory of the driv- 
ing system and let y 6 M d denote the trajectory of the 
response system. Hence, the full phase space that de- 
scribes the evolution of the two systems is M 2d . If z 
denotes a trajectory in the full phase space, then the 
synchronization manifold for identical systems is defined 
byW s = {z= (x,y) |x = y}. 

Next, assume that the driving trajectory evolves on 
an m < d dimensional invariant manifold, Wd- This 
manifold may, or may not, be stable to perturbations. 
For example, if the driving trajectory is an unstable fixed 
point then m = 0, and Wd is unstable. If the driving 
trajectory is an unstable limit cycle then m = 1, and 
Wd is unstable. However, if the driving trajectory is a 
dense orbit on a chaotic attractor then m < d and Wd 
is stable. (For this case we are not saying that Wd is 
the attractor. Rather, Wd is the manifold that contains 
the attractor.) For all of these cases it is possible for Wd 
to occupy zero volume in R d . Finally, if the system is 
Hamiltonian and the driving trajectory is chaotic then 
m = d and Wd occupies positive volume in R d . 

In addition, the dynamical system itself may have an 
invariant manifold, W° C IR d (which we call an invariant 
manifold of the system.) An invariant manifold of the 
system may contain an infinite number of trajectories. 
For example, the Ott-Sommerer system discussed above 
has an invariant manifold defined by y — Vy — 0, and 
a chaotic attractor exists in this manifold |2l|| . (This 
chaotic set is attracting for points within the manifold 
defined by y = v y = 0.) Therefore, it is possible to 
consider a case where Wd is restricted to an invariant 
manifold of the system, Wd Q W°. However, in general 
this is not the case. An important point to keep in mind is 
that the driving trajectory defines Wd, and each different 
driving trajectory (in general) defines a different Wd- 

For identical systems, the response system has an iden- 
tical m dimensional manifold, Wr. The manifolds Wd 
and Wr are not the same because Wd lives in the phase 
space of the driving system while Wr lives in the phase 
space of the response system. 

In the full phase space of the combined system these 
manifolds can be defined by W D = {z = (x, y) |x 6 Wd}, 
and W* R = {z = (x,y) | y S Wr}. (See Fig. [l| for a 
schematic picture of these manifolds.) These manifolds 
intersect in the full phase space. The intersection is de- 
fined by W D C\W* R = {z = (x, y)|x e W D and y e Wr}, 
(see the dark solid curve in Fig. 12). Since this defini- 
tion does not require x = y it is clear that, although 
W D p| W R contains part of the synchronization manifold, 
W* D f| W* R / W s . Furthermore, W s contains all possible 
driving trajectories. Therefore, if Wd is restricted to an 
invariant manifold of the system (Wd W°) then some 
trajectories within W s are not in W D . 

Therefore, when we discuss the linear stability of syn- 
chronization for a given driving trajectory we are ac- 
tually discussing the stability of the manifold given by 




FIG. 12. Various manifolds and intersections which arise in 
the phase space associated with two identical systems which 
are coupled together. 



W = W D nmn^ s = {* = (x,y) xeW^aiidye 
W R and x = y} to perturbations that are transverse to 
this manifold. The manifold associated with this inter- 



section is labeled p in Fig 12, Although it is shown as 



a dot one must keep in mind that the dimension of the 
manifold, W, ranges from a low of zero to a high of d. If 
W D P| W r P| W s is stable for all possible driving trajec- 
tories we say that W s is stable. 



APPENDIX C: GEOMETRICAL 
INTERPRETATION OF THE STABILITY 
CRITERIA 

In this appendix we discuss a geometrical interpreta- 
tion of the stability criteria discussed in Sections [lj and 
III. Equations (10) and ( |l2| ) have geometrical interpre- 
tations which can be used to design couplings that yield 
stable synchronous motion. One can think of the ele- 
ments of DE(0) as living in a d 2 dimensional parameter 
space. The right hand side of Eq. ( |l0| ) defines a function 
in this parameter space. Furthermore, this function can 
be used to define a family of surfaces in this parameter 
space. Explicitly, the family of surfaces £b is defined 
by (HP -1 [B(x;i)]P||) = C B , where C B is a constant. 
In a similar fashion the left has side defines a different 
function and family of surfaces. Explicitly, the family of 
surfaces Ea is defined by — 5ft[Ai] = Ca, where Ca is a 
constant. 

Now consider a fixed driving trajectory. For this tra- 



17 



C„=3 




FIG. 13. Various surfaces Ea and Eb- The dashed line is 
the boundary of the region of parameter space that yields lin- 
early stable synchronization. Any location in parameter space 
"above" this line will yield linearly stable synchronization. 



jectory, each value of Cb (CO corresponds to a partic- 
ular surface from the family Eb (^a). The boundary of 
the region of parameter space that yields linearly stable 
synchronization is given by the intersections of surfaces 
from Ea with surfaces from Eb where Cb = Ca- These 
intersections form a surface in the parameter space. Any 
intersection associated with Ca > Cb will result in lin- 
early stable synchronization onto the driving trajectory. 
Furthermore, all of these intersections will reside on one 
side of the surface defined by the boundary (see Fig. 13). 
Because Ca and Cb are (in general) arbitrary real num- 
bers, and the reals are dense, we know that the region 
associated with linearly stable synchronization will typi- 
cally occupy positive volume in the phase space. 

An immediate consequence of this approach is that if 
none of the surfaces from Ea intersect surfaces from Eb 
then the rigorous criterion for synchronization can not 
be met. 

An explicit examples of this analysis is the Ott- 
Sommerer example discussed above. For this example 
Eqs. (0) and (^8[) are equivalent to the stability crite- 
rion of Eq. ( pp| ) . Begin with the diagonal driving example 
and consider a fixed driving trajectory. For this exam- 
ple the parameter space is 1R and C (from Eq. (|T§|)) is a 
constant, independent of e. Because C is constant, Cb 
can have only one value and there is only one surface in 
the family Eb- Furthermore, the fact that the right right 
hand side of Eq. (jl^) is independent of e implies that Eb , 
is the entire parameter space. In contrast Ca can take 
on any value. The left hand side of Eq. JL7| ) indicates 
that a particular value of Ca corresponding to a particu- 
lar value of e. Therefore, Ea is a family of points in the 
parameter space ]R. 

For these surfaces the boundary of the region of param- 
eter space that yields linearly stable synchronization (the 
intersection of Ea and Eb associated with Ca = Cb ) is 



a point in JR. Intersections associated with Ca > Cb 
are also points. Finally, the region associated with lin- 
early stable synchronization (the intersection of all sur- 
faces from Ea with surfaces from Eb such that Ca > Cb) 
is an interval which occupies positive volume in M (see 
Eq.©). 

The final example is driving with position, and the po- 
sition is permitted to drive both the position and velocity 
equations. The parameter space for this example is M 2 . 
Now, assume the driving trajectory is fixed. The deriva- 
tion that lead to Eq. (|l|) shows that, Eb is a family of 
parabolas in M 2 . Selecting a value for Cb is equivalent 
to selecting a value for C, which is equivalent to selecting 
a particular parabola from this family. 

The left hand side of Eq. @ indicates that Ea is a 
family of lines in R 2 . Selecting a value for Ca selects a 
particular line from this family. Thus, if one has chosen a 
particular value for C then the boundary of the parame- 
ter space region that yields stable synchronization is the 
intersection of a line and a parabola. For this example 
this intersection always exists. Intersections of lines as- 
sociated with Ca > Cb (for fixed C) sweep out a portion 
of the parabola. Therefore, for fixed Cb (equivalently, 
fixed C) the portion of the parabola that is beyond the 
line associated with Ca = Cb corresponds to values of 
the coupling parameters that produce stable synchronous 
motion. 

Now recall that the value of C is arbitrary, and differ- 
ent values of C yield different parabolas. For each dif- 
ferent value of C the rigorous criterion chooses a portion 
of a parabola. Therefore, the region of parameter space 
that satisfies the rigorous criterion is the one swept out 
by portions of the parabolas for all possible values of C. 
This region occupies positive volume in the parameter 
space, M 2 . 



APPENDIX D: RELATIONSHIP TO CONTROL 
THEORY 

In this appendix we discuss relationship between our 
results and those found in control theory, and alternate 
formulations of the rigorous criterion. The relationship 
between synchronization and control has been know for 
some time [ ^8|38[ |. In these papers synchronization is 
typically discussed as an example of feedback control (see 
the review by Chen and Dong ||). 

In basic control theory one is often interested in the 
following equation 



dw 
~dt 



A(t)w + B(t)u(t) 



where u(t) is a time dependent scalar input (the drive). 
In this equation we are invited to think of A as the un- 
controlled dynamics and B(t) is a vector which couples 
the input to the uncontrolled dynamics. An important 
question in control theory is, can one find inputs, w(i), so 
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that the w(t) = in finite time. The answer to this ques- 
tion is often determined by constructing a controllability 
matrix 

G(t,t ) = f \](t,s)B(s)B t (s)V' f (t,s)ds, (Dl) 

J t 

where U(t,to) = expfj^ A(s)ds}. It is known that if 
G(t,to) has full rank (det(G) ^ 0) then one can find 
inputs so that w(t) = |3jj |. 

We are interested in the following linearized equation 
of motion 



dw 
~dt 



[DF(f) + DE(0)] w. 



For the coupling we have considered the inputs are feed- 
backs, u(t) = C a w a where C is a constant vec- 
tor (which typically has only one component.) There- 
fore, a controllability approach would involve determin- 
ing a B so that G(t, to) has full rank and then using 
DE(O) a0 = B a Cf3 to obtain DE(0). 

The difficulty with the controllability approach is that 
\J(t,to) has a complicated time dependence. Avoiding 
this complication is one of the motivations for considering 
the equation 



dw 
~dt 



= [A + B(t)]w, 



(D2) 



where A and B are defined by Eqs. (||and (^). In this for- 
mulation of the problem the coupling strengths (DE(0)) 
are part of A, we think of A as the uncontrolled dynam- 
ics, and B(t)w are driving inputs. In effect the rigorous 
criterion of Eq. (10) says that if the eigenvalues of A are 
placed sufficiently far into the right half plane (as a re- 
sult of feedback coupling) from the fluctuation that result 
from the chaotic driving will not be sufficient to destabi- 
lize the fixed point at w = 0. In this sense our criterion 
is similar in spirit to pole placement control theory. 

We remark that the controllability matrix can not be 
immediately calculated for Eq. (D2) because B is a func- 
tion that maps input space into state space while B is a 
function that maps state space into itself. In Eq. (D2) 



the function that plays the role of B is part of the matrix 
A. 

We believe that a different interpretation of synchro- 
nization (still from control theory) could be quite promis- 
ing. Walker and Mees interpret synchronization as an 
example state estimation and the observer problem [p| . 
This view of synchronization says that synchronization 
permits one to estimate the complete state of a system 
when the full state vector is not observable. In practice 
one could imagine coupling the output from a physical 
system to a model of the system using coupling that is 
guaranteed to result in stable synchronization. (Here we 
rely on the fact that synchronization is robust to model- 
ing errors and noise The full state of the physical is 



then the same as the state of the model. For linear sys- 
tem the observer problem and the control problem can 
can be shown to be duals of each other EJuJ . Thus, 
many of the results proved for observability of a dynami- 
cal system are corollaries of theorems proved for control. 
In this paper a rigorous and easy to calculate criterion 
for synchronization is given. 



APPENDIX E: ANOTHER RIGOROUS 
CRITERION 

In this appendix we describe a derivation that leads 
to a rigorous criterion similar to the one derived in Sec- 
tion [n| If one does not transform to coordinates given 
by the eigenvectors of A then Eq. (||) becomes 

w(t) = U(f,to)w(to)+ / V(t ) s)B(s)w(s)ds, 

where U(t, to) — exp[A(t — to)]. This equation leads to 
Eq. (0) as a sufficient condition for linear stability of syn- 
chronization with K replaced by B. It is straightforward 
to shown if U(i, to) = exp[A(t — to)] then a rigorous suf- 
ficient condition for linear stability of synchronization is 



K[A 1 ]>a(||B|| 



where 



C = 



1/3 



1/2 



(El) 



(E2) 



For this criterion the terms involving the coupling 
strength and the time average have been decoupled. This 
could greatly simplify the effort required for some calcu- 
lations. 

As an example of this consider the Lorenz system. The 
eigenvectors of A are 



e± 



N± 



[(a-e 4 ), (a + ei+A±), 0], 



where N± are the following normalizations 



N± = [{a - 



E4 



|<7 + ei + A : 



,2-1 V2 



Also, B(x) is given by 



B(x) 





— z + (z) —x + (%) 

y - (y) x ~ (x) o 



Clearly, calculating ||P~ 1 B(x)P|| is a complicated pro- 
cedure. However, the rigorous condition of Eq. (El) can 
be easily satisfied. It is possible to show that if A± are 
complex then 
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c 



(r-(z)-e 3 + l) : 



+ ex - 1 - e 2 ] 2 - 4(r - (z) - e 3 )(a - e 4 ) 



1/2 



Now consider off diagonal y-driving (ei =£3 = 0) and 
fixed C. To satisfy the rigorous condition of Eq. (El) 
while keeping C fixed one needs 



(e a + 1 - <t) 



4(r-(z» 
e 2 > -((j + l) + 4C. 



(^-(^)-l) 2 
4C 2 (r- (z)) 



|B||), 



where C is arbitrary. The second equation is a parabola 
which exists for all values of £2, and all driving trajecto- 
ries. The first equation is a line which is guaranteed to 
intersect the parabola. Therefore the rigorous condition 
of Eq. (pi]) can be satisfied for this type of driving. 
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Measure Type 


Drive Type 


Approximate 


Test 




Numerical Test 






£ (C) 




£ (c) 








(c) 






mm 




c max 




mm 




c max 




X 


0.1975 




4.998 




0.1975 




4.998 


Fixed Point 


V 


0.1976 




oo 




0.1976 




oo 




z 


9.0 




zzo 




9.0 
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X 


0.04 




4.8 




0.54 




3.7 


Period 1 


y 


0.04 




oo 




0.36 




OO 




z 


F 




r 






None 






X 


0.07 




4.9 




0.32 




4.3 


Period 2 


y 


0.07 




oo 




0.29 




oo 




z 


F 




F 






None 






X 


0.11 




4.9 




0.20 




4.5 


SBR 


y 


0.11 




oo 




0.18 




oo 




z 


F 




F 






None 




TABLE I. Results of numerical tests on the Rossler system. 


In this table F implies 


that this type of driving fails the test, 


while none implies 


that synchronization did not occur. 


















Rossler System 








Lorenz System 




Measure Type 


(x) 


(y) 


<*> 




» 


(y) 




<^> 


Period 1 


0.2770 


-1.385 


1.385 












54.81 


Period 2 


0.2246 


-1.123 


1.123 












54.94 


SBR 


0.1649 


-0.8245 


0.824f 


> 










54.82 


TABLE II. Numerical estimates for averages along driving trajectories. 


For the Rossler system 


(x) = —a 


(y) and (z) = - (y). 






Lorenz System 












Measure Type 


Drive Type 


Approximate 


Test 




Numerical Test 










(c) 




e (c) 




(c) 






mm 




c max 




mm 




c max 




X 


590 




OO 




590 




oo 


Fixed Point 


y 


59 




OO 




59 




oo 




z 


F 




F 






None 






X 


42 




oo 




17 




oo 


Period 1 


y 


4.2 




oo 




4.1 




oo 




z 


F 




F 






None 






X 


42 




oo 




18 




oo 


Period 2 


y 


4.2 




oo 




4.0 




oo 




z 


F 




F 






None 






X 


42 




oo 




19 




oo 


SBR 


y 


4.2 




oo 




4.0 




oo 




z 


F 




F 






None 




TABLE III. Results of numerical tests 


on the Lorenz system. 


In this table F implies 


that this type of driving fails the test, 


while none implies 


that synchronization did not occur. 


















Ott-Sommerer System 










Measure Type 


<^> 


<<? (ll) ) 




<l5 (±) -<3 ( 


X) )D 






-<. 9 (ll) )l) 


Period 1 


-1.223 


-6.307 




0.4769 






5.142 


Period 2 


-3 


-7.767 




1.678 








10.30 


SBR 


Q 

-o 


-7.038 




1.714 








7.856 


TABLE IV. 
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Drive Measure 


Drive Type 


Ott-Sommerer System 
Rigorous Tests Approximate Test 

-L II -L II 


Numerical Test 
J- II 


Period 1 


Diagonal All 


0.812 


4.04 








0.01 


0.28 




Diagonal only x 


F 


F 








0.01 


0.62 




Diagonal only v 


F 


F 








0.01 


0.69 




Off Diagonal x 


0.95 


0.95 








0.002 


0.15 









99.6 








0.04 


0.55 


Period 2 


Diagonal All 


1.89 


7.34 








0.13 


0.99 




Diagonal only x 


F 


F 








0.25 


2.6 




Diagonal only v 


F 


F 








0.24 


2.8 




Off Diagonal x 


0.95 


0.95 








0.05 


0.5 






8.47 


416 








0.9 


1.75 


SBR 


Diagonal All 


1.92 


5.87 








0.05 


0.52 




Diagonal only x 


F 


F 








0.09 


1.9 




Diagonal only v 


F 


F 








0.10 


1.5 




Off Diagonal x 


0.95 


0.95 








0.05 


0.5 






8.95 


193 








0.7 


1.75 



TABLE V. Results of numerical tests on the Ott-Sommerer system. In this table F implies that this type of driving fails the 
test. For off diagonal driving the first number listed is either ei or £2 while the number listed below it is either £3 or £4. 
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